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The energy spectrum of spin-orbit coupled states of individual sub-surface boron acceptor dopants 
in silicon have been investigated using scanning tunneling spectroscopy (STS) at cryogenic temper¬ 
atures. The spatially resolved tunnel spectra show two resonances which we ascribe to the heavy- 
and light-hole Kramers doublets. This type of broken degeneracy has recently been argued to be 
advantageous for the lifetime of acceptor-based qubits [Phys. Rev. B 88 064308 (2013)]. The 
depth dependent energy splitting between the heavy- and light-hole Kramers doublets is consistent 
with tight binding calculations, and is in excess of 1 meV for all acceptors within the experimen¬ 
tally accessible depth range (< 2 nm from the surface). These results will aid the development of 
tunable acceptor-based qubits in silicon with long coherence times and the possibility for electrical 
manipulation. 


Dopant atoms in silicon are attractive candidates for 
spin-based quantum computation. Recent studies have 
demonstrated long coherence-times for both ensembles 
of bulk donors [T] and individual donors [2] in silicon. 
Meanwhile, rapid progress in scanning tunnelling mi¬ 
croscopy (STM) based lithography has paved the way 
towards atomically precise placement of dopant atoms 
[3]. While phosphorous donors in silicon remains among 
the most compelling candidates for dopant-based quan¬ 
tum computing to date, other impurity systems have re¬ 
cently drawn considerable attention. In particular, boron 
acceptors could provide a pathway towards electrically 
addressable spin-qubits via spin-orbit coupling [4] ana¬ 
logues to electrically driven spin manipulation in gate 
defined electron 016 ] and hole [ 7 ] quantum dots in III-V 
materials. Compared with other spin-orbit qubits, ac¬ 
ceptors in silicon have several advantages: gates are not 
required for hole confinement and each qubit experiences 
the same confinement potential; furthermore the hole- 
spin decoherence, due to the nuclear spin bath, can be 
effectively eliminated by isotope purification of the silicon 
host. 

Unlike the ground state of donor-bound electrons in 
silicon, the acceptor-bound hole ground state is four-fold 
degenerate, reflecting the heavy-hole/light-hole degener¬ 
acy of the silicon valence band. Recent theoretical work 
has suggested the regime of long lifetimes for acceptors 
with a four-fold degenerate ground state is only accessi¬ 
ble for small magnetic fields [4]. Interestingly, this work 
also suggests that symmetry breaking due to strain or 
electric fields could yield longer-lived qubits based on 
acceptor-bound holes at higher magnetic fields. The sym¬ 
metry breaking perturbation of biaxial strain [4] renders 
the lowest two (qubit) levels within a Kramers degener¬ 
ate pair, such that they do not directly couple to elec¬ 
tric fields. Quantum confinement could provide a similar 


form of protection. In this Letter we demonstrate that 
the symmetry-reduction of a potential boundary renders 
the lowest two levels Kramers degenerate and heavy-hole 
like. Recent transport spectroscopy studies of an indi¬ 
vidual acceptor embedded in nano-scale transistors have 
shown that for acceptors ^10 nm away from an interface 
the bulk-like four-fold degeneracy is maintained [8] . Here 
we demonstrate that the presence of a nearby interface 
(<2 nm) lifts the four-fold degeneracy of the acceptor- 
bound hole ground state by investigating the energy spec¬ 
trum and wavefunctions of individual sub-surface boron 
acceptors using scanning tunnelling spectroscopy (STS). 

Recent low-temperature STM/STS studies have suc¬ 
cessfully probed the energy spectrum and wavefunctions 
of individual impurity atoms in Si EHni and GaAs 
nans]. These studies have revealed the influence of the 
semiconductor/vacuum interface on the ionization en¬ 
ergy of sub-surface dopants laiisi, the spatially resolved 
structure of the dopant wavefunction noHia and the 
mechanisms for charge-transport through these dopants 
El HB E]. Here, we use STS to measure the energy 
difference between the heavy-hole and light-hole states 
of individual B acceptors less than 2 nm away from the 
surface. The sample is prepared by repeated flash an¬ 
nealing a 8 X 10^^ cm“^ boron doped Si(lOO) substrate 
in ultra-high vacuum. The flash annealing not only pro¬ 
vides an atomically flat surface but also yields a region of 
low doping at the interface due to out-diffusion of the B 
impurities. As a result, the near-interface acceptors are 
weakly coupled to the bulk acceptors, which form a va¬ 
lence impurity band (VIB). Together with the STM tip, 
the sub-surface acceptors can be considered as a double¬ 
barrier system for single-hole tunneling E] (see Fig. 1(a)). 
As illustrated in the inset of Figure 1(b), individual ac¬ 
ceptors are identified as protrusions at Vb = —1.5 V due 
to enhancement of the valence band density of states El • 
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(a) Tip Vacuum Sample 



FIG. 1. (a) Schematic energy diagram of the tunneling pro¬ 
cess. Applying a positive sample bias (eVb) results in an up¬ 
ward tip-induced band bending (TIBB). When the energy of a 
state of an isolated acceptor near the interface is pulled above 
the Fermi level (Ef) there is a peak in the conductance due 
to resonant tunneling from the valence impurity band (VIB) 
through the acceptor state into the tip. In this schematic de¬ 
scription, the top layer of silicon is depleted of dopants. The 
measured acceptors are well-isolated from each other, as well 
as the underlying nominally 8 x 10^® cm“^ doped region, de¬ 
noted by VIB. (b) Current and (c) differential conductance 
measured directly above the isolated acceptor (open circles) 
as a function of bias voltage. The differential conductance 
within the band gap is fitted to the sum (solid line) of two 
thermally broadened Lorentzian line-shapes corresponding to 
resonant tunnelling through the lowest two acceptor states. 


Figure l(b,c) shows the current and differential con¬ 
ductance measured over a single sub-surface acceptor. 
For the bias voltage < 0 V the charge transport is 
dominated by holes tunnelling directly from the Si va¬ 
lence band to the tip. Similarly, for V 5 > 1.1 V trans¬ 
port is dominated by holes tunnelling from the tip to 
the Si conduction band. The observed features in the 
band-gap, i.e. for 0 V< V 5 < 1.1 V, can be attributed 
to holes tunnelling from the VIB through the localized 
acceptor states to the tip. Upon increasing the bias volt¬ 
age from 0 V to 1.1 V the energy levels of the subsurface 
acceptor states shift due to tip induced band bending 
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FIG. 2. (a) Spatially resolved differential conductance map 
as a function of tip position and bias voltage measured over a 
sub-surface acceptor (different from Fig. 1). (b,c) Probability 
density map of the lowest two acceptor states. 


(TIBB). Each time an acceptor level enters the bias win¬ 
dow, which is defined by the Fermi level of the tip and the 
Si substrate, this opens an additional channel for trans¬ 
port resulting in a stepwise increase in the current (Fig. 
1 (b)) and a corresponding peak in the differential con¬ 
ductance as shown in Fig. 1(c). The conductance peaks 
are only observed in the presence of an acceptor, and 
for each measured acceptor we observe two conductance 
peaks in the band-gap which we attribute to the lowest 
two acceptor states entering the bias window. 

A spatially resolved differential conductance map is 
shown in Fig. 2(a). The amplitude of the differential 
conductance measured over the acceptor is a measure 
of the probability density of the impurity wavefunction 
(Fig. 2(b,c)). The probability density maps of the low¬ 
est two acceptor states have the same spatial extent and 
symmetry, i.e. neither state shows anti-nodes in the sur¬ 
face plane. The similarity between the probability den¬ 
sity map of the two conductance peaks is consistent with 
what is expected for ±3/2 and ±1/2 states. Both are ex¬ 
pected to have the observed s-like envelopes, with similar 
spatial extents m- In contrast, a two-hole state would 
have a larger spatial extent. When two states have the 
same envelope wavefunctions the energy difference be¬ 
tween these two states can only arise from the difference 
in pseudo-spin, i.e. from non-degenerate heavy-hole and 
light-hole states. In the remainder of this Letter we will 
discuss how the perturbation of the subsurface acceptor 
wavefunctions by an interface lifts the degeneracy be¬ 
tween the heavy- and light-hole Kramers doublets that 
make up the four-fold degenerate ground state of an un¬ 
perturbed (bulk) acceptor. 

Based on the tunneling spectra of the single acceptor 
we can now extract the energy splittings of the states, 
and show that they are in reasonable agreement with 
tight binding predictions for the s-like 3/2 and 1/2 states. 
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and too small to be associated with charging processes 
where a second hole is bound to the acceptor. If the 
lever arm a = dEi/edVi that couples the chemical po¬ 
tential to the bias voltage is known, the voltage Vi for 
which a localized state i is brought into resonance is a 
direct measure for the eigenenergy Ei of this state. We 
determine a by fitting the conductance to a thermally 
broadened Lorentzian [16] (see Fig. 1(b)): 
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where is the Boltzmann factor, T the temperature 
(we assume T = 4.2 K for the fit), Vi the voltage cor¬ 
responding to the center of the conductance peak, h 
Planck’s constant and F the sum of the tunnel-in and 
tunnel-out rates. The interface induced splitting is given 
by 6E = a{V±i /2 - 143 / 2 )- 

The distance of the subsurface acceptors to the inter¬ 
face is measured from the spectral shift of the valence 
band edge |9|: 
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where d is the acceptor depth and Q = e(e^ — esi)/{ey + 
esi) is the image charge due to the mismatch between 
the dielectric constants €y and esi of the vacuum and sil¬ 
icon, respectively. The onset voltage for tunneling from 
the valence band Vy is determined by finding the volt¬ 
age axis intercept of the linear extrapolation of the nor¬ 
malised conductance curve at its maximum slope point 
dl- The spectral shift of the valence-band edge is fitted 
to equation using the depth d of individual acceptors 
and the modified dielectric constant Q as two, indepen¬ 
dent, fitting parameters. The obtained values for Q for 
all measured acceptors agree within experimental error 
with the expected value Q = e(e^ — esi)/{^v + ^Si) fol¬ 
lowing the classical half-space approach and experimen¬ 
tal values that have previously been reported from STM 
experiments nansj. 

Figure 3 shows the measured energy splitting 5E be¬ 
tween the lowest two acceptor states as a function of ac¬ 
ceptor depth for six different acceptors demarcated by the 
solid diamonds. These results are compared with a fully 
atomistic tight binding (TB) simulation marked by the 
grey circles. The tight-binding Hamiltonian of 1.4 million 
silicon atoms with a boron acceptor was represented with 
a 20-orbital sp3d5s* basis per atom including nearest- 
neighbor and spin-orbit interactions. An acceptor was 
represented by a Coulomb potential of a negative charge 
screened by the dielectric constant of Si and subjected 
to an onsite cutoff potential Uq [20]. The model provides 
an accurate solution for the single-hole eigenstates of a 
bulk acceptor, and an acceptor near an interface. The 
magnetic field is represented by a vector potential in a 
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FIG. 3. Experimental and theoretical depth dependence of 
the interface induced energy splitting. The vertical error bars 
on the energy-level splitting are proportional to the confidence 
of the fit of a. Errors in Vi are negligible. The horizontal 
error bars are proportional to the confidence of the fit of d 
which decreases for acceptors further from the interface as the 
spectral shift of the valence-band maximum becomes smaller. 


symmetric gauge and entered through a Peierls substitu¬ 
tion. The full TB Hamiltonian is solved in NEM03D [21] 
by a parallel Block Lanczos algorithm, and the relevant 
low energy acceptor wavefunctions are obtained. From 
the calculated magnetic field dependence of the accep¬ 
tor states (not shown) we infer that the lowest Kramers 
doublet has a predominant heavy-hole (±3/2) character 
whereas the second Kramers doublet has a predominant 
light-hole (±1/2) character. 

From the combination of the spatial and spectral mea¬ 
surements, we can conclude that the observed states 
belong to the s-like manifold which contains a ±3/2 
Kramers doublet and a ±1/2 Kramers doublet. Not only 
are both of the envelopes s-like, but the few meV energy 
scale of the splitting is too small to be associated with 
higher excited states (at 25 meV) [22] or charging tran¬ 
sitions (47 meV) [23] of an acceptors. We note that the 
charging energy near the surface could fall as low as 20 
meV m, but this is still too large for our results. 

We observe an increase of SE for acceptors closer to 
the interface, as predicted by the TB calculations. Im¬ 
portantly, all acceptors studied showed an energy split¬ 
ting in excess of 1 meV. For comparison, the electric field 
required to obtain a splitting of 0.5 meV would be 40 
MV/m [24] , much greater than the electric field required 
for field ionization 5 MV/m [25]. Our results therefore 
demonstrate that the presence of an interface provides 
an effective way to energetically isolate a single Kramers 
doublet that could serve as the working levels of a spin- 
qubit. Such a qubit could benefit from elimination of the 
nuclear spin bath by isotope purification, and could pro¬ 
vide a new route towards an electrically controllable spin 
qubit. 
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